Is the misalignment of the Local Group velocity and the 2MRS dipole typical in a 

ACDM & halo model? 
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We predict the acceleration of the Local Group generated by the 2MASS Redshift Survey (2MRS) 
within the framework of ACDM and the halo model of galaxies. We show that as the galaxy 
fluctuations derived from the halo model have more power on small scales compared with the mass 
fluctuations, the misalignment angle between the CMB velocity vector and the 2MRS dipole is 
in reasonable agreement with the observed 21°. This statistical analysis suggests that it is not 
necessary to invoke a hypothetical nearby galaxy or a distant cluster to explain this misalignment. 

PACS numbers: 98.56 



I. INTRODUCTION 

In the late 1980s and early 1990s, the near alignment 
of the IRAS and the optical dipoles with the Local 
Group velocity relative to the cosmic microwave back- 
ground (CMB, [U, i, i, i, i, i, 0, i) were used to verify 
the gravitational instability picture of structure forma- 
tion and the scale of convergence was used to assess the 
validity of the then popular Standard Cold Dark Matter 
(SCDM) model (0 hereafter (LKH) and [3, hereafter 
JVW). 

The distance at which most of the Local Group (LG) 
peculiar velocity is generated {the convergence depth) has 
remained a contentious issue for the past twenty five 
years. The observed misalignment between the vectors 
derived from the galaxy survey and the CMB has been 
attributed to several factors: 

1. The linear perturbation theory of density fluctua- 
tions is correct only to first order density contrast, 
5. There may be contributions to the LG dipole 
from small scales which would cause gravity and 
the velocity vectors to misalign. 

2. The selection effects of the surveys from which 
these vectors are calculated will increase the shot 
noise-error especially at large distances causing 
misalignments. 

3. There may be uncertainties in the assumptions in 
galaxy formation and clustering. For example, the 
mass-to-light ratios might differ according to type 
and/or vary with luminosity or the galaxy/cluster 
biasing might be non- linear and/or scale depen- 
dent. 

4. There may be a significant contribution to the LG 
dipole from structure further away than the maxi- 
mum distance of the surveys studied. 

5. The direction of the LG dipole may be affected by 
nearby galaxies or distant clusters behind the Zone 
of Avoidance (ZoA), which are not sampled by most 
surveys. 

The estimates of the convergence depth are highly sensi- 
tive to the type of data and the analysis performed. For 



example, the convergence depth estimated from galaxy 
clusters consistently disagree with the estimates from 
galaxy surveys (e.g. see and the references therein). 
Moreover, the treatment of the ZoA plays an impor- 
tant factor in dipole determinations. Exclusion of bright 
galaxies changes the dipole misalignment direction con- 
siderably (from 21° to 14°, [ll|). These bright galaxies 
lie within the ZoA and are excluded from most dipole 
analyses. This also highlights the vital role non-linear 
dynamics introduced by nearby objects play in dipole 
considerations. The analysis of the convergence of the 
dipole is further complicated by the redshift distortions 
on small and large scales which introduce systematic er- 
rors [13. 

Recently, the new dipole determinations from the 
2MASS Redshift Survey (hereafter 2MRS, [O, [H, [H, [3 
brought up once again the question of alignment and con- 
vergence. For example, [lB| attempted to fix the observed 
misalignment by placing a hypothetical object behind 
the ZoA and [131 tried to mitigate the non-linear and 
shot noise effects by removing nearby objects in order to 
increase the correlation between the dipole determined 
from the CMB and the 2MRS. Furthermore, two recent 
studies of the peculiar velocity surveys indicated a large 
convergence depth. ^18^ found a strong bulk flow on scales 
out to 300 /i-^Mpc u sing peculiar velocity data from an 
X-ray cluster sample. [19| reported a large-scale bulk flow 
beyond 50 h~^Mpc from a comprehensive compilation of 
peculiar velocity data. Both these papers claimed that 
these convergence depths, determined by their studies, 
were difficult to explain within the framework of stan- 
dard ACDM model of cosmology. 

Here we take a different approach, updating the LKH 
and JVW calculations to include two major developments 
in cosmology: the ACDM concordance model (replacing 
the old SCDM and a halo model for the galaxy power 
spectrum (replacing the naive linear biasing model). In 
the next section, we consider the expected misalignment 
between the peculiar velocity of the LG and the accelera- 
tion vector estimates. We expect that the ACDM model 
will give a poorer misalignment compared to the SCDM 
model since it has a larger coherence length. Adding the 
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halo model for galaxy biasing will further increase the 
power of density fluctuations on small scales and there- 
fore the misalignment angle. In Section 3, we consider 
the convergence of the acceleration from the 2MRS sur- 
vey. In Section 4, we show that the observed convergence 
is highly probable in our model. In the last section we 
discuss limitations of our calculations and future work. 



II. ALIGNMENT 

The most popular mechanism for the formation of 
large-scale structure and motions in the Universe is the 
gravitational growth of primordial density perturbations. 
According to this paradigm the peculiar acceleration vec- 
tor g(r) induced by the matter distribution around posi- 
tion r is related to the mass by 



oo 

g(r) = GpJ dh'S^ir')^^^ 



(1) 



where p is the mean matter density and dm (r) — {pm (r) — 
p)/p is the density contrast of the mass perturbations. 
In linear theory, the peculiar velocity field, v(r), and the 
peculiar acceleration vector, g(r), should be parallel with 
a constant of proportionality. 

For v(r), we use a value derived from the CMB dipole. 
Using the first year of data from WMAP, [23| and the 
revised values of the motion of the Sun relative to the LG 
a velocity derived by 21], we find a LG velocity relative 
to the CMB of vlg = 627± 22 kmsec~\ towards {Ilg = 
273° ± 3°,6lg = 29° ± 3°). The LG dipole, d(r) (in 
kmsec"^), used as a proxy for g(r), is estimated from a 
galaxy survey and can be written as a summation of all 
the galaxies in that survey (Equation 13 in fllj): 



d(r) = -Si > WL,Sir, , 

PlOl ^ 



(2) 



where ^(V) is the luminosity function and ium.i is the 
minimum luminosity at a distance (see [TT| for a full 
derivation). The Poisson shot noise of d is estimated as 



N 



(4) 



For a whole sky survey, the non-zero elements of the 
covariance matrix are given as follows: 



Cdd = ^(d.d) 



and 



Cvv= 3(v.v) 



and 
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Wl{k)Pm„,{k)dk^ (6) 
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where Pgg and Pmm 



Wv(fc)VKd(fc)Pg,„(fc)dfc3, 
(7) 

denote the galaxy and the dark mat- 
ter power spectrum, respectively and Pgm denotes the 
galaxy-dark matter cross spectrum. 

The expected misalignment angle between v and d is 
(JVW and LKH): 



p{9\v)^sm{e) exp (-1/(202)) 



exp (2/V2) 
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1 + erf (y/V2) 



(8) 



where 9 is the angle between the two vectors, 
9^ = [Vl — (? I c\l (y I cV^) (in radians) and c = 
[C|^/(CvvC'dd)]^^2 and y = cos{9)/9^. In the case where 
0* is small, the probability distribution p(cos(0)|v) tends 
to bi-variate Gaussian and 9^ measures the expected scat- 
ter. 



where pi, is the luminosity density, &i is the luminos- 
ity bias factor introduced to account for the dark matter 
haloes not fully represented by the galaxies (hereafter we 
assume 6^ — 1). The index 7 appears in the growth 
function / = d\n6/d\na « [22] and is an important 
diagnostic for distinguishing between ACDM and alter- 
native models of gravity. It has little dependence on the 
cosmological constant !23'l, and a slight dependence on 
the dark energy equation of state parameters w 2J| . Re- 
cent refined calculations predict 7 = 0.55 for the ACDM 
concordance model, and 7 = 0.69 [1^ for a particular 
modified gravity model, DGP braneworld gravity [26| . 

In equation [51 Si is the flux of galaxy i and w^. is the 
luminosity weighting function written as: 



r.L^{L)dL' 



(3) 



A. Window Functions 

We assume the LG is a sphere with radius tlg = 
l.Qh~'^ Mpc and the LG velocity is produced by the mat- 
ter distribution out to horizon. In Fourier space, the win- 
dow function of the velocity vector v, determined from 
the CMB dipole, varies as 1/k: 



W^r{k) = 



jo(fc'"LG) 



(9) 



For the flux- weighted dipole W^{k) can be written as: 



Wjj(fc) = 47r / r^drW^{r)ii{kr) 

_ j0(fc?'LG) _ j0(fc?-„ax) 

~ k k ' 



(10) 
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where Wjj(r) = and rmax is the maximum radius of 
the survey which we adopt as r^ax = 130h~^ Mpc. We 
want to compare our findings with the observed 2MRS 
values and beyond this distance the shot-noise effects in 
the 2MRS becomes too large. 



B. Calculating the Power Spectra 

The halo model for non- linear clustering gained popu- 
larity during the past few years [13, [H, [2§| . This model 
assumes that most matter in the Universe is structured 
as virialised halos and correlation function on large scales 
depends on the clustering of these halos whereas non- 
linear clustering of objects are related to the internal 
structure of the individual halos. This implies that the 
non-linear power spectrum P„i{k) can be divided into a 
quasi-linear, inter halo term P2h{k) and a non-linear in- 
tra halo term Pih{k): 



Pnl{k) = P2h{k) + Plhik). (11) 



C. The Dark Matter Power Spectrum 

Within the halo model framework, the dark matter 
power spectrum can be written as (e.g. f30|). 

^in^(fc) - /dM n(M) (^ly |u(k|M)p, (12) 

and 

H'^ik) = P^\lik) y dM n(M)b(M) (^j^ u(k|M) 

^ as) 

where the integrals are over the halo mass, M and p is 
the background density of the universe. We use transfer 
functions of [31] to calculate P^™(fc). In the expressions 
above, n{M) is the halo mass function describes the num- 
ber density of haloes of mass M, b(M) is the halo bias 
factor and u{k\M) is the Fourier transform of the halo 
profile, p{r). 

For n(M) and b(M), we assume the forms proposed by 



n(M)dM = ^i{v)Av, 



M 



1/2 



-qi./2 



(14) 



(15) 



where q — 0.707, p — 0.3, and the normalisation 
Af, K, 0.3222. The mass variable is defined as v = 
((5sc(z)/ct(M))^, where 5sc{z) is the hnear-theory predic- 
tion for the present day overdensity of a region under- 
going spherical collapse at redshift z ((5sc(0) ~ 1.68) and 
a{M) is the r.m.s. variance of the present day linear 



power spectrum in a spherical top-hat which contains an 
average mass M . 

The degree of bias is also a function of the halo mass: 



6(M) = 1 + 



qv -I 2p/Ssciz) 



5,,iz) l + iqjy)P' 
We also use the [33|] dark matter halo density profile 



(16) 



p{r) 



Ps 



{r/r,){l+r/r,y 



{r < Tvir) 



(17) 



where rg is the characteristic scale radius and ps provides 
the normalization. The profile is truncated at the virial 
radius r^ii, which is obtained from the halo mass. We 
parametrize the profile in terms of the concentration pa- 
rameter c = r^iy/rg. The normalization for the mass M 
is 



M 



p{r)4:Trr'^dr = Anpgr^ 



In (1 + c) 



1 



(18) 

We assume that the concentration parameter c depends 
on halo mass M and redshift in a manner calibrated by 
numerical simulations ( [33. [35l| ) : 



c(M,z = 0) = 11 I — 



-0.13 



(19) 



where Mq is obtained from equation by setting the vari- 
ance of the linear power spectrum (t(Mo, 0) = i5sc = 1.68. 
For our adopted cosmological model we obtain Mq = 

101213 /j-IMq. 



D. The Galaxy Power Spectrum 

The components of the galaxy power spectrum are sim- 
ilar to the dark matter power spectrum given above 

Pir(fc)= /dMn(M)»i_l»|,(k|M)|2, (20) 



"gal 



and 



PS^\k) = P'J^Lik) 



dM n(M)b(M) 



(N|M) 

hgal 



u(k|M) 



(21) 

Here, {N\M) and (N{N - 1)|M) are the first and second 
factorial moments of the halo occupation distribution, 
P{N\M), respectively, rigai is the average number density 
of galaxies. The halo occupation distribution (HOD) is a 
parametrized description of how galaxies populate dark 
matter haloes as a function of the halo mass M. We 
adopt simple models motivated by results from simula- 
tions and semi-analytic calculations (e.g. [3^ ) 



{N\M) 



{M/Mo 




i/3 



(M > Mcut) 
(M < Mcut) 



(22) 
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where Mq and /3 are free parameters and Mcut is deter- 
mined by matching figai to the value derived from the 
2MRS (ngai = 7.7 X IQ-'^h^Mpc^^). The second moment 
of P{N\M) is: 

^ (23) 

where a = log^JM/h-^10^^MQ. 

For our halo model parameters, we use Mq = 
IO^ i^/i-^Mq and (3 = 0.85, derived for 2MASS galaxies 
(see Table 1 of [131 )• These values are quite typical for 
red galaxies and varying the parameters slightly does not 
affect the results. 

We note that on large scales the galaxy power spectrum 
takes the simple form: 

Pgai(A;)=62pr.(fc) + l/^igai (24) 

The last term is the shot noise term, which we found to 
give the very r.m.s. contribution to the shot-noise in the 
dipole as the empirical estimate (Equation [4]) . 

E. The Galaxy-Dark Matter Cross Power 
Spectrum 
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FIG. 1: The probability distribution functions for misalign- 
ment angle G between the galaxy survey dipole and the CMB 
velocity vectors. The dashed-dotted line denotes the misalign- 
ment angle calculated using P(k)|j5„ only for SCDM model 
((1) in Table 1), the dashed line is calculated using ACDM 
and P(k)nJjn ((3) in Table 1). The solid line is calculated us- 
ing Equations [HE] and [7]((4) in Table 1). The dotted line is 
the same as the solid line but using P(k)gg = P(k)gg ((5) in 
Table 1). 



Using the expressions above, the galaxy dark matter 
cross power spectrum components can be written as 



PiiJ^)(fc)= /dMn(M)Mm|t,(k|M)P, 
P ngai 



(25) 



and 



dM n(M)b(M)^^i^u(k|M) 



Ilgal 



dMn(M)b(M) ( — ] u(k|M)] 



M 



(26) 



III. CONVERGENCE 

We now turn to the question of the convergence of the 
acceleration vector derived from a galaxy survey given 
the velocity of the CMB. We recalculate the conditional 
probability,p(d|v) given in equation [8] within a series of 
successively larger concentric spheres centred on the LG. 
In this case, the value of the covariance matrix element 
Cvv remains the same and we recalculate C^v and Cdd 
using a series of window functions: 



io(fcfLG) joikR) 



(27) 



where R is the outer limit of each shell. 

The expected dipole velocity given the CMB dipole is 
then written as: 



d(i?) 



-VCMB 



(28) 



and the '1 sigma' scatter in the velocity is |d(i?)| ± votr^ 
with = C, 



dd 



IV. RESULTS 

In Figure [U we plot a family of p{9\v) curves given by 
equation [51 The dashed-dotted line is calculated using 
the old SCDM model with = 1.0, as = 0.6, the Hub- 
ble constant Hq = 50kmsec~^Mpc~^ and the spectral in- 
dex n = 1. For all other power spectra, we use WMAP5 
fiducial parameters fim = 0.26, fb = 0.17, erg = 0.79, 
Ho = 72kmsec-^Mpc~^ and n = 1 

We first consider the standard linear and non-linear 
ACDM dark matter power spectrum models with a lin- 
ear biasing model of 6 = 1. We obtain almost identi- 
cal curves (dashed line) with a scatter of 9^, =7°. Both 
curves peak at around 8° . Changing the radius of the LG 
to tlg = 5h~^ Mpc has also a very small effect on the re- 
sults. The narrowest scatter in the probability curve and 
the smallest expected misalignment angle is obtained us- 
ing the SCDM model {9^ = 6°, dashed-dotted hue). This 
is expected since SCDM has a smaller coherence length 
than the ACDM model. 

We then calculate the power spectra using the halo 
model for galaxies as well as for dark matter. The power 
spectrum shown by the black solid line includes both 1- 
halo and the 2-halo terms. The scatter for this model 
is the largest 6** = 14° and the probability curve peaks 
at the largest misalignment angle of the plotted curves 
(13.5°) . In the galaxy halo model, the galaxies are more 
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FIG. 2: The linear dark matter (Piin, solid line), halo model 
dark matter (Pdm, dashed line), galaxy-dark matter (Pgm, 
dashed-dotted line) and galaxy-galaxy (-Pgai, dotted line) 
power spectra. The power spectra are normalized to be equal 
on very large scales (very small k). 

clustered than the dark matter (see Figure [2|) and thus 
the covariance matrix elements are much larger than the 
ones we obtain with the linear bias models. The halo 
model is used to describe all haloes within a survey. Thus 
the 2-halo term represents the clustering of all haloes 
on large scales and the 1-halo term is related to the in- 
ternal clustering of each halo in the survey. We also 
want to model the dipole induced by clusters of galax- 
ies rather than galaxies themselves. For this, we plot the 
probability function calculated using 2-halo term only 
{Pgg{k) = P2h{k), the dotted line) which peaks around 
10° and gives 61* = 11°. When we add the 1-halo term 
(solid line), we are adding the galaxy component to our 
model. This increases the shot-noise which causes the 
CMB and the galaxy survey dipoles to decorrelate as con- 
firmed by our value for 9^. Changing the radius of the 
LG to tlg = 5h^^ Mpc mitigates the non-linear effects, 
reducing the scatter, 6*, = 10°. 

We also calculate the r.m.s. velocity of our LG models, 
< Vlg >^/^= (3C™)i/2. For our ACDM cosmology with 
as = 0.79, we obtain 412 and 430 kmsec^^for the linear 
and non-linear power spectra, respectively. We note the 
observed LG velocity amplitude, 621 kmsec"^, is larger 
than these r.m.s. values. 

We summarize the results presented in Figure [T] in Ta- 
ble 1. We present the probability obtaining the actual 
observed angle (21°) and the integrated probability for 
9 to lie below the observed value which provides a fre- 
quentist test. The two columns indicate that the ob- 
served dipole is highly unlikely in both the linear and 
non-linear models which do not include a non-linear bias 
model for galaxies. Neither are we able to approximate 
the observed value without including the 1-halo term to 
the galaxy power spectrum. But, using the halo model 
for both galaxies and the dark matter, results in a mean 
close to that of 21°. 

In Figure [31 we plot the expected dipole velocity given 
the CMB dipole (Equation [551 solid line) and its scatter 
(dashed line). In the top panel, we use only the hnear 
power spectrum (the first model in Table 1) and the bot- 
tom panel is the velocity convergence with the full halo 
model treatment (the fourth model in Table 1). The 



plots in the left panel are the predicted amplitudes of 
the acceleration on the LG as a function of distance. The 
right plots show the convergence of the angle between the 
galaxy survey dipole and the CMB dipole. We also plot 
the observed values of the acceleration amplitude and the 
misalignment angle from the 2MRS survey obtained us- 
ing the fiux- weighted selection function ([11]) as an illus- 
trative example. There are many determinations of the 
dipole misalignment ang le from the 2MRS data ranging 
from 14° Jdj] to 50° [lj| depending on the type of anal- 
ysis performed. The near-infrared fiux weighted dipoles 
are more robust that the number weighted schemes be- 
cause they closely approximate a mass-weighted dipole, 
bypassing the effects of redshift distortions and require no 
preferred reference frame. Our analysis holds for either 
scheme or the choice of reconstruction technique since we 
account for the correct data vector for selection effects 
and compare like with like. 

The 2MRS acceleration grows rapidly out to « 
50h^^ Mpc after that it fiattens off. The misalignment 
angle drops to 12° at « 5Qh~^ Mpc but then increases 
at larger distances. This behaviour is not seen for either 
of the models. This is not unexpected, however, since 
the increase in amplitude at « 5Qh~^ Mpc and than its 
consequent decrease at « 60h~^ Mpc is attributed to the 
tug-of-war between the Great Attractor and the Pisces- 
Perseus Supercluster which are not accounted in the sta- 
tistical (r.m.s.) models. 

In theory, one can equate the velocity inferred from the 
CMB measurements with the value derived from a galaxy 
survey and obtain a value for Q,^. In practice, however, 
the galaxy surveys do not measure the true total veloc- 
ity due their finite depth and the predicted amplitude 
of the reconstructed velocity depends the combination of 
matter density and the biasing parameter b which is usu- 
ally not known to great accuracy. However, one would 
expect fim and b inferred from the LG velocity to have 
observationally viable values. Figure [3] shows that the 
amplitudes for both of the models are similar to observa- 
tions, but the non- linear model has higher amplitude due 
to higher small-scale power and is closer to the observed 
values. In both models, more than half of the predicted 
LG velocity signal is generated within « 30h~^ Mpc. The 
non-linear galaxy model decorrelates the angles and in- 
corporates shot-noise effects and closely resembles the 
observed convergence. 



V. DISCUSSION 

We have compared estimates for of the acceleration as 
derived from a 2MRS like survey and the motion of the 
LG in a ACDM universe both within and without the 
halo model framework. We showed that a combination 
of linear theory with the halo model predicts convergence 
and misalignment very much in agreement with the ob- 
served 2MRS dipole. 

We should note our analysis assumes linear perturba- 
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TABLE I: Dipoles constrained by the velocity of the local group, the characteristic misalignment angle ^* is given in degrees. 
The observed misalignment angle Oobs ~ 21° is calculated from the 2MRS survey (Erdogdu et al. 2006b). The assumed 
cosmological parameters {a8,ho,^m,fb) are (0.79,0.72,0.26,0.17) for ACDM and (0.6,0.5,1.0,0.044) for SCDM. 



Model 


e,{deg) P(6'obs)/max 


[P{e\v)) Percent below 6'obs 


(1) SCDM, P^'^(fc) 


6.3 


0.02 


99.6 


(2) ACDM, P^^ik) 


7.2 


0.08 


98.1 


(3) ACDM, Pdl(fe) 


7.2 


0.08 


98.1 


(4) ACDM, Halo model for Galaxies 13.5 


0.78 


67.2 


(5) ACDM, P^^{k)^P^^ik) 


10.4 


0.45 


85.3 



? from LG (Mpc/h) 



' from LG (Mpc/h) 
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60 00 100 120 
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FIG. 3: Top Panel: The predicted LG velocity within a se- 
ries of successively larger concentric spheres centred on the 
LG given by Equation 1281 and the 1-sigma scatter. Bottom 
Panel: The misalignment angle as a function of distance. 
Left: The convergence is calculated using linear dark matter 
power spectrum only ((2) in Table 1). Right: The conver- 
gence is calculated using non-linear galaxy power spectrum 
and galaxy-dark matter cross spectrum (Case (4) in Table 4) . 



tion theory but we are using a non-linear galaxy power 
spectrum, so our treatment is incomplete. A more consis- 
tent way would be to use N-body simulations (e.g. [3§|) 
and semi-analytic models galaxy formation but this is be- 
yond the scope of the work presented in this paper. Here, 
we tested if a non-linear power spectrum model with a 
more sophisticated biasing formula could predict the ob- 
served acceleration better than a linear power spectrum. 
Indeed, our results indicate the flow in our location in the 
universe is quite typical for a ACDM model. A non-linear 
velocity-acceleration relation will increase the misalign- 
ment angle Recently, Lavaux et al. used the 



2MRS reconstructed peculiar velocity field to estimate a 
LG velocity misalignment of « 50° ± 22° (at 95% confi- 
dence). They suggest the dipole directions should agree 
to within 25° at 95% confidence. From Figure 3 (right 
panel) , the expected misalignment angle is ~ 29° ± 7° at 
130 Mpc. Our analysis should be viewed as a con- 
servative estimate of the expected misalignment error. 
Thus, even a dipole misalignment is as high as reported 
by jlj] is reconcilable with a ACDM model given their 
estimated error. 

We would like to emphasize that we have assumed 
that the accelerations are estimated from a whole sky 
survey. It is possible to incorporate the effect of ZoA 
(see Appendix of LKH) which would increase the scatter 
substantially. We note that we are able to reproduce 
the observed values very well without modeling miss- 
ing structures in the masked regions as suggested by 
We further note that many dedicated searches of 
the ZoA did not detect such massive structures (see e.g. 

[u, El, m, m, m, m, e m, m, iio, im ) . our analysis 

should be understood in the r.m.s. sense. It does not 
pick out directions. As noted by others (e.g. [3), most 
of the LG misalignment lies along the direction of the 
galactic centre, the region of most obscuration, and this 
is another route for exploring the misalignment. 

To summarize, we presented a fresh analysis of the 
fundamental problem of the origin of motion on the Lo- 
cal Group. We conclude that the newly observed 2MRS 
dipole is in accord with what we would expect in a ACDM 
universe within the context of the halo model. 



VI. ACKNOWLEDGEMENTS 

We thank Michal Chodorowski, Gert Hutsi and Renee 
Kraan-Kortewcg for helpful discussions. PE would like to 
thank UCL for its hospitality. OL acknowledges the sup- 
port of a Wolfson Royal Society Research Merit Award. 



[1] Yahil A., Walker D. & Rowan-Robinson M., 1986, ApJ, [2] Meiksin A. & Davis M., 1986, AJ, 91, 191 
301, LI 



7 



[3] Harmon R.T., Lahav O. & Meurs E.J.A., 1987, MNRAS, 
228, 5 

[4] Villumsen J.V. & Strauss M.A., 1987, ApJ, 322, 37 

[5] Lahav O., Rowan-Robinson M. & Lynden-Bell, 1988, 

MNRAS, 234, 677 
[6] Lynden-Bell D., Lahav O. & Burstein D., 1989, MNRAS, 

241,325 

[7] Strauss M. A., Yahil A., Davis M., Huchra J. P., Fisher 

K., ApJ, 1992, 397, 395 
[8] Webster M., Lahav O. & Fisher K., 1997, MNRAS, 

287,425 

[9] Lahav O., Kaiser N., Hoffman Y., 1990, ApJ, 352, 448 
[10] Juszkiewicz R., Vittorio N., Wyse R., 1990, ApJ, 352, 
408 

[11] Erdogdu P., et al., 2006a, MNRAS, 368, 1515 

[12] Kaiser N., 1987, MNRAS, 227, 1 

[13] Erdogdu P., et al., 2006b, MNRAS, 373, 45 

[14] Lavaux G., Tully R. B., Mohayaee R., Colombi, S., 2008, 

arXiv:astro-ph/0810.3658v2 
[15] Crook A., Silvestri A., Zukin P., 2009, arXiv:astro- 

ph/0906.2411 
[16] Loeb A., Narayan R., 2008, MNRAS, 386, 2221 
[17] Chodorowski M. J., Coiffard J., BiUcki M., Colombi S., 

Ciecielag P., 2008, astro-ph/0706.0619 
[18] Kashlinsky A., Atrio-Barandela F., Kocevski D., Ebeling 

H., 2008, ApJ, 686, 49 
[19] Watkins R., Feldman H. A., Hudson M. J., 2009, MN- 
RAS, 392, 743 
[20] Bennett C.L., et al. , 2003, ApJS, 148, 1 
[21] Courteau S. & Van Den Bergh S., 1999, AJ, 118, 337 
[22] Peebles P. J. E., The Large- Scale Structure of the Uni- 
verse, Princeton: Princeton Univ. Press, Princeton, N.J., 
USA 

[23] Lahav O., Lilje P. B., Primack J. R., Rees M. J., 1991, 

MNRAS, 251,' 128 
[24] Wang L., Steinhardt P. J., 1998, ApJ, 508, 483 
[25] Linder E. V., Cahn R. N., 2007, APh, 28, 481 
[26] Dvah C, Gabadadze G., Porrati M., 2000, PhLB, 485, 

208 

[27] Ma C, Fry J. N., 2000, ApJ, 543, 503 

[28] Peacock J.A., Smith R.E., 2000, MNRAS, 318, 1144 



[29] Seljak U., 2000, MNRAS, 318, 203 
[30] Cooray A., Sheth R. K., 2002, Phys. Rep., 371, 1 
[31] Eisenstein D. J., Hu W., 1998, ApJ, 496, 605 
[32] Sheth R. K., Tormen G., 1999, MNRAS, 308, 119 
[33] Navarro J. F., Frenk C. S., White S. D. M., 1997, ApJ, 
490, 493 

[34] Bullock J. S., Kolatt T. S., Sigad Y., SomerviUe R. S., 
Kravtsov A. V., Klypin A. A., Primack J. R., Dekel A., 
2001, MNRAS, 321, 559 

[35] Zehavi I. et al., 2004, ApJ, 608, 16 

[36] Kravtsov A. V., Berlind A. A., Wechsler R. H., Klypin 
A. A., Gottlober S., AUgood B., Primack J. R., 2004, 
ApJ, 609, 35 

[37] Lin Y., Mohr J.J., Stanford S.A., 2004, ApJ, 610, 745 
[38] Dunkley J. et al. 2009, ApJ Supp., 180, 306 
[39] Davis M., Strauss M.A., Yahil A., 1990, ApJ, 372, 394 
[40] Davis M., Nusser A., Willick J., 1996, ApJ, 473, 22 
[41] Ebeling, H., MuUis, C. R., Tully, R. B. 2002, ApJ, 580, 
774 

[42] Ebeling H., Kocevski D., Tully R. B., MuUis C. R. 2005, 

Nearby Large-Scale Structures and the ZoA, 329, 83 
[43] Kocevski, D. D., Ebeling, H., MuUis, C. R., 
Tully, R. B. 2005, ArXiv Astrophysics e-prints, 
arXiv:astro-ph/0512321 
[44] Kocevski D. D., Ebehng H., 2006, ApJ, 645, 1043 
[45] Roman, A. T., Takeuchi, T. T., Nakanishi, K., Saito, M. 

1998, PASJ, 50, 47 
[46] Wakamatsu K., Malkan, M. A., Nishida M. T., Parker 
Q. A., Saunders W., Watson F. G. 2005, Nearby Large- 
Scale Structures and the ZoA, 329, 189 
[47] Hasegawa T., et al. 2000, MNRAS, 316, 326 
[48] Meyer M. J., et al. 2004, MNRAS, 350, 1195 
[49] Henning P. A., Kraan-Korteweg R. C, Stavely-Smith L. 
2005, Nearby Large-Scale Structures and the ZoA, 329, 
199 

[50] Kraan-Korteweg R. C, Lahav 2000 A&A Rev, 10, 211 
[51] Kraan-Korteweg R. C, Shafi N., Koribalski B., Staveley- 

Smith L., Buckland P., Henning, P. A., Fairall, A. P. 

2007, larXiv:0710.1795l 



